OPTIMAL LONG TERM INVESTMENT MODEL 
WITH MEMORY 
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Abstract. We consider a financial market model driven by an R. n -valued 
Gaussian process with stationary increments which is different from Brownian 
motion. This driving noise process consists of n independent components, and 
each component has memory described by two parameters. For this market 
model, we explicitly solve optimal investment problems. These include (i) 
Merton's portfolio optimization problem; (ii) the maximization of growth rate 
of expected utility of wealth over the infinite horizon; (iii) the maximization 
of the large deviation probability that the wealth grows at a higher rate than 
a given benchmark. The estimation of paremeters is also considered. 



1. Introduction 

In this paper we study optimal investment problems for a financial market model 
with memory. This market model M. consists of n risky and one riskless assets. 
The price of the riskless asset is denoted by Sq (t) and that of the iih risky asset 
by Si(t). We put S(t) = (Si(t), . . . , S n (t))', where A' denotes the transpose of 
a matrix A. The dynamics of the R™-valued process S(t) are described by the 
stochastic differential equation 



r x — ^ n 

dSi(t) = Si(t) m(t)dt + =1 ^j(t)dY 3 (t) 

(1.1J 



t > 0, 

Si(0) = Si, i = l,...,n, 

while those of So(t) by the ordinary differential equation 

(1.2) dS<,(t) = r(t)S (t)dt, t>0, S (0) = 1, 

where the coefficients r(t) > 0, fJ-i(t), and <Jij{t) are continuous deterministic func- 
tions on [0, oo) and the initial prices Sj are positive constants. We assume that the 
n x n volatility matrix a(t) = (<Tjj(t))i<jj< n is nonsingular for t > 0. 

The major feature of the model Ai is the R"-valued driving noise process Y(t) — 
(Yi(t), . . . ,Y n (t))' which has memory. We define the jth component Yj(i) by the 
autoregressive type equation 

where W(t) = (W\{t), . . . , W n (t))' , t e R, is an R"-valued standard Brownian 
motion defined on a complete probability space (O, J 7 , P) , the derivatives dYj(t)/dt 
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and dWj (t) / dt are in the random distribution sense, and pj 's and qj 's are constants 
such that 



(1.4) < qj < oo, —qj < p < oo, j = 1, . . . , n 

(cf. Anh and Inoue 1 ). Equivalently, we may define Yj(t) by the moving-average 
type representation 



(1.5) Y j (t) = W J (t) 



Pj 



e -(l:i+P^ S - U UWj(u) 



ds, t <E R 



(see Examples 2.12 and 2.14]). The components Yj(t), j — 1, . . . , n, are Gaussian 
processes with stationary increments that are independent of each other. Each Yj (t) 
has short memory that is described by the two parameters pj and qj . In the special 
case pj = 0, Yj(t) reduces to the Brownian motion Wj(t). Driving noise processes 
with short or long memory of this kind are considered in £Q, Anh et al. and 
Inoue et al. [201, for the case n = 1. 
We define 

T t := o{a{Y{s) : < s < t) UjV) , t > 0, 

where N is the P-null subsets of T. This filtration (J-t)t>o is the underlying 
information structure of the market model Ai. From l|1.5|l . we can easily show that 
(Y(t))t>o is a semimartingale with respect to (Tt) (cf. Section 3]). In particular, 
we can interpret the stochastic differential equation l|l.l|) in the usual sense. In 
actual calculations, however, we need explicit semimartingale representations of 
Y(t). It should be noticed that l|1.5fl is not a semimartingale representation of Y(t) 
(except in the special case Pj = 0) . For, Wj (t) involves the information of Y 3 (s) with 
s < and vice versa. The following two kinds of semimartingale representations of 
Y(t) are obtained in (21 Example 5.3] and |2()l Theorem 2.1], respectively: 

ft r 



(1.6) Y j (t)=B j (t) 



o L^o 
* r 



(1.7) 



Y j (t)=B j (t) 



kj(s, u)dYj(u) 



lj(s, u)dBj(u) 



ds, t > 0, j = 1, . . . , n, 



ds, t > 0, j = 1, . . . , n, 



where, for j = 1, . . . ,n, (Bj(t)) t >o is the so-called innovation process, i.e., an R- 
valued standard Brownian motion such that 

a{Y 3 {s):0<s<t) = a(Bj (s) : < s < t), t> 0. 

Notice that Bj's, are independent of each other. The point of (|1.6|) and (|1.7|) is that 
the deterministic kernels kj(t, s) and lj(t, s) are given explicitly by 

{2q 3 +p 3 )eii s -p 3 e-i^ 



(1.8) 

(1.9) 
with 

(1.10) 



kj (t, s) = Pj (2 qj +Pj) (2g . +p . )2e ^_ p g e - git , 
lj(t, s) = e~^ +q i ){t - s) l 3 {s), 0<s<t, 



0<s<t, 



lj(s) :=p 3 



2p 3 q 3 



(2qj +p 3 ) 2 e 2 <ii s -p) 



s > 0. 



We have the equalities 



(1.11) 



k 3 {t,s)dYj(s) 



lj(t,s)dB 3 (s), t>0, j = l, 



Many authors consider financial market models in which the standard driving 
noise, that is, Brownian motion, is replaced by a different one, such as fractional 
Brownian motion, so that the model can capture memory effect. To name some 
related contributions, let us mention here Comte and Renault [SJ , Rogers [SD] , 
Heyde QB], Willinger et al. |32), Barndorff-Nielsen and Shephard 0, Barndorff- 
Nielsen et al. g], Hu and 0ksendal [IS). Hu et al. [TSJ, Elliott and van der Hoek 0, 
and Heyde and Leonenko |1 7| . In most of these references, driving noise processes 
are assumed to have stationary increments since this is a natural requirement of 
simplicity. Among such models, the above model M. driven by Y(t) which is a 
Gaussian process with stationary increments is possibly the simplest one. One 
advantage of M. is that, by the semimartingale representations and l|1.7|l of 
Y(t), it admits explicit calculations in problems such as those considered in this 
paper. Another advantageous feature of the model A4 is that, assuming (t) = , 
real constants, we can easily estimate the characteristic parameters pj, qj and cr^ 
from stock price data. We consider this parameter estimation in Appendix C. 

For the market model M., we consider an agent who has initial endowment 
x G (0, oo ) and invests fi(t)X I,,r (i) dollars in the ith risky asset for i = 1, . . . ,n 
and [1— X^iLi 7T i(t)]X x,7T (t) dollars in the riskless asset at each time t, where X x ' w (t) 
denotes the agent's wealth at time t. The wealth process X x '*(t) is governed by 
the stochastic differential equation 



(1.12) 



dSo(t) , v^n us dSi(t) 



X*,T(t) L ^i=i n >\ S (t) ^i=i n ; Si(t) 

Here, we choose the self-financing strategy ir(t) = (7Ti(i), . . . ,7r„(i))' from the ad- 
missible class 

{7T is an R™-valued, progressively measurable^ 
= {n(t))o<t< T : . T ,^ n2 , M \ 

process satisfying J Q \\Tr(t)\\ at < oo a.s. J 

for the finite time horizon of length T 6 (0, oo), where || • || denotes the Euclidean 
norm of R". If the time horizon is infinite, we choose n(t) from the class 

A := {(7r(t)) t >o : (7r(t)) <t<T S A T for every T e (0, oo)} . 

Let a S (— oo, 1) \ {0} and c 6 R. In this paper, we consider the following three 
optimal investment problems for the model M.: 

(PI) V(T,a) := sup -E[(X^(T)) a ], 

tt£A t a 

(P2) J(a) := sup lwasup-^- log E[(X^(T)) a } , 

TteA T-hx al 

(P3) I{c) := sup limsup - logP [X s '* (T) > e cT ] . 

The goal of Problem PI is to maximize the expected utility of wealth at the 
end of the finite horizon. This classical optimal investment problem dates back 
to Merton We refer to Karatzas and Shreve (23 and references therein for 

work on this and related problems. In Hu et al. |T2], this problem is solved for 
a Black-Scholes type model driven by fractional Brownian motion. In Section [2 
assuming pj > for j — 1 , . . . , n, we explicitly solve this problem for the model 
M. . Our approach is based on a Cameron-Martin type formula which we prove in 
Appendix A. This formula holds under the assumption that a relevant Riccati type 



equation has a solution, and the key step of our arguments is to show the existence 
of such a solution (Xemma l2.1[l . 

The aim of Problem P2 is to maximize the growth rate of expected utility of 
wealth over the infinite horizon. This problem is studied by Bielecki and Pliska 
[H], and subsequently by other authors under various settings, including Fleming 
and Sheu HUES, Kuroda and Nagai [221, Pham [23123, Nagai and Peng [2"Y| . 
Hata and Iida ^Sj , and Hata and Sekine |14l ITS] . In Section we solve Problem 
P2 for the model A4 by verifying that a candidate of optimal strategy suggested 
by the solution to Problem PI is actually optimal. In so doing, existence results 
on solutions to Riccati type equations (Lemmas 12.11 and 13.5(1 play a key role as 
in Problem PI. The result of Nagai and Peng [27] on the asymptotic behavior of 
solutions to Riccati equations, which we review in Appendix B, is also an essential 
ingredient in our arguments. 

The purpose of Problem P3 is to maximize the large deviation probability that 
the wealth grows at a higher rate than the given benchmark c. This problem is 
studied by Pham |28l I2H| . in which a significant result, that is, a duality relation 
between Problems P2 and P3, is established. Subsequently, this problem is studied 
by Hata and Iida ^2] and Hata and Sekine |14l ITS] under different settings. In 
Section 21 we solve Problem P3 for the market model M. In the approach of 
|28U29| . one needs an explicit expression of J (a). Since our solution to Problem P2 
is explicit, we can solve Problem P3 for M. using this approach. As in .2 8", I55| . our 
solution to Problem 3 is given in the form of a sequence of nearly optimal strategies. 
For c < c with certain constant c, an optimal strategy, rather than such a nearly 
optimal sequence, is obtained by ergodic arguments. 



2. Optimal investment over the finite horizon 

In this section, we consider the finite horizon optimization problem PI for the 
market model M.. Throughout this section, we assume a 6 (— oo, 1) \ {0} and 

(2.1) < qj < oo, < pj < oo, j = 1, ... ,n. 

Thus pj > rather than pj > —qj (see Remark 12.61 below 1 . 

Let Y(t) = (Yi(i),...,Y n (t))' and B(t) = {B x {t), . . . , B n {t))' be the driving 
noise and innovation processes, respectively, described in Section^ We define an 
R n - valued deterministic function X(t) = (Ai (<),... , X n (t))' by 

(2.2) \(t):=a-\t)[n(t)-r(t)l], t > 0, 

where 1 := (1, . . . , 1)' G R". For the kernels kj(t, s)'s in i|1.8fl . we put 

k(t, s) :— diag(fci(t, s), . . . , k n {t, s)), < s <t. 
We denote by = (£i (t),..., £ n (t))' the R"-valued process f Q k{t, s)dY(s), i.e., 

(2.3) &(<):= / k j (t,s)dY j (s), t>0, j = l,...,n. 

Jo 

By <|1.1[1 . (|l-2p. I|1.6fl . and l|1.12|l . the wealth process X x ' % (t) evolves according to 
dX x ^(t) 



= r{t)dt + Tr'{t)a{t) [X(t) - f (t)] dt + 7r'(t)a(t)dB(t), t > 0, 



whence, by the Ito formula, we have, for t > 0, 



X x ' n (t) = a; exp 



(2.4) 



r(s) + n'(s)a(s) (A(,s) - £(,)) - -\\a'(s)Tr(s) 



ds 



-!- / 7r'(s)o-(s)dB(s) 
/o 



We define an R- valued process Z(t) by 
Z{t) := exp 



£ > 0. 



{A(s) -£(*)}' dB(a)-- / |]A( S )-e( S )fd S 

Since A(i)— is a continuous Gaussian process, the process is a P-martingale 
(see, e.g., Example 3(a) in Liptser and Shiryayev |23l Section 6.2]). We define the 
R-valued process (F(t))o<t<T by 

T(t) := E [Z^(T)\ Tt] , < t < T, 

where (3 is the conjugate exponent of a, i.e., 

(1/a) + (1//3) - 1. 

Notice that < < 1 (resp. -oo < < 0) if — oo < a < (resp. < a < 1). In 
view of Theorem 7.6 in Karatzas and Shreve Chapter 3], to solve Problem PI, 
we only have to derive a stochastic integral representation for T(t). 
We define an R-valued P-martingale K(t) by 

rt pa rt 

K(t) := exp 



-13 ^ {\(8)-t(8)YdB(a) -Lj \\ X {s) -H{s)\\ 2 ds 



t > 0. 



Then, by Bayes' rule, we have 



r(t) = e 



K (T) exp i --0(1-0) ^ \\\(8)-t(a)\\'d8 



K(t)E 



rap j --0(1-/3) ^ ||A( S )-e( S )|| 2 ^ 



for i G [0, T], where E stands for the expectation with respect to the probability 
measure P on (Q, T T ) such that dP/dP = K(T). Thus 

' T ||A(,s)fd S 
I - J i 

(2.5) 



T(t) = Z?(t) exp | -1/3(1 -0)jf 



x E 



1 



exp 



0(1-/3) / (U(s)\\ 2 ~2X'(s)as))ds 



We are to apply Theorem lA.ll in Appendix A to (|2.5() . By Ijl.lip . the dynamics 
of are described by the n-dimensional stochastic differential equation 

(2.6) d£(t)=-(p + q)£(t)dt + l(t)dB(t), t>0, 

where p := diag(pi, . . . ,p n ), q := diag(gi, . . .,q n ), and l(t) := diag(Zi(i), . . .,l n (t)) 
with L(tYs as in fl~TU|l . Write := B(t)+0 J*[X(s) — ^(s)]ds for t G [0,T]. Then 
P(i) is an R"-valued standard Brownian motion under P. By (|2.6|l . the process 
evolves according to 

(2.7) d£(t) = \p(t) + b(t)£(t)]dt + l(t)dB(t), t > 0, 



where p(t) = (pi(t), . . .,p n (t))', b(t) = diag(&i(t), . . . ,b n (ij) with 

(2.8) Pj (t) :=-plj it)\j(t), t>0, j = l,...,n, 

(2.9) 6i(*):=-(Pi+*)+/3ii(*), *>0, j = l,...,n. 

By Theorcm lA.ll in Appendix A, we are led to consider the following one-dimensional 
backward Riccati equations: for j = 1, . . . , n 

(2.10) Rj(t) - $(t)i$(i) + 2bj{t)Rj(t) + 0(1-0)=O, < t < T, Rj (T) = 0. 
The following lemma, especially (iii), is crucial in our arguments. 

Lemma 2.1. Let j G {1, . . . , n}. 

(i) //pj = 7 t/ien H2.10(l /ias a unique solution Rj(t) = Rj(t;T). 

(ii) // — oo < a < 0, i/ien <|2.10ll /ias a unique nonnegative solution Rj(t) = 
Rjit-T). 

(iii) //pj > and < a < 1, then (|2.10() has a unique solution Rj(t) = Rj(t; T) 
such that Rj(t) > bj(i)/l?(t) for t £ [0, T\. 

Proof, (i) If pj — 0, then 12.10J1 is linear, whence it has a unique solution. 

(ii) If — oo < a < 0, then 0(1 — 0) > 0, so that, by the well-known result on 
Riccati equations (see, e.g., Fleming and Rishel ^3 Theorem 5.2] and Liptser and 
Shiryayev |231 Theorem 10.2]), l|2.10|l has a unique nonnegative solution. 

(iii) When pj > and < a < 1, write 

(2.11) oi(t) := if (t), a 2 (t);=bj(t), a 3 := 0(1-0), t > 0. 
Then the equation for P(t) := Rj(t) — [a 2 (t) / a\(t)\ becomes 

(2.12) P(t) - ai {t)P 2 {t) + a A (t) = 0, 0<£<T, 
where 

a 2 (t) + ai(t)a 3 d \ a 2 (t) 



ai{t) - 0l (t) ' dt 

Since dlj(t)/dt > and < 0, we see that 



ai (i) 



02(t) 



2(Pj + gj )-%(*) ^j/, Wn 



_<H(t)_ 

We write a 2 (t) + ai(t)a 3 as 

(1 - 0) [fa + q 3 f - {( Pj + q,) - i,(<)} 2 ] + [fa + qj) - Ijit)} 2 , 

which is positive since < lj(t) < pj. Thus (t) > 0, so that 12.12|l has a unique 
nonnegative solution P(t) = P(t;T). The desired solution to 12.1011 is given by 
R j (t) = P(t) + [a 2 (t)/a 1 (t)}. □ 

In what follows, we write Rj(t) = Rj(t;T) for the unique solution to l|2.10|l in 
the sense of Lemma l2.1l Then R(t) := diag(i?i(t), . . . , R n (t)) satisfies the backward 
matrix Riccati equation 

R(t) - R(t)l 2 (t)R(t) + b(t)R(t) + R(t)b(t) + 0(1 - 0)I n = 0, < t < T, 
< ~ 2 ' 13 ' 1 R(T) = 0, 



(2.14) 



(2.15) 



where /„ denotes the n x n unit matrix. For j = 1, . . . ,n, let Vj(t) = Vj(t;T) be 
the solution to the following one-dimensional linear equation: 

+ hit) - (*;!>,-(*) + [3(1 - (3)Xj(t) - Rj(t;T) Pj (t) = 0, 

0<i<T, Vj(T) = 0. 

Then w(i) = v(t;T) :— (vi(t;T), . . . ,v n (t;T))' satisfies the matrix equation 

v(t) + [b(t) - l 2 (t)R(t; T)]v(t) + (3(1 - f3)X{t) - R(t; T)p(t) = 0, 

< t < T, v(T) = 0. 

We put, for j = 1, . . . , n and (t, T) e A, 

(2.16) g 3 {t;T) := v^t; T)l](t) + 2 Pj (t)v {t:T) - l^R^T) - (3{1 - f3)X 2 j (t), 
where 

(2.17) A := {(t, T) : < T < oo, < i < T}. 
We are now ready to give the desired representation for T(t). 

Proposition 2.2. Write 

(2.18) ip(t) :=T(t)[-(3X(t) + {(3-l{t)R(t;T)}t(t) + l(t)v(t;T)}, < t < T. 
Then, for t G [0,T], we /iaue T(t) = T(0) + f*ip'(s)dB(s) with 

1 r T n 

Proo/. It follows from l(23|l. (|2~7|) . fTIS) . (|2~15|) and Theorem [Q that 



(2.19) 



r(0) = exp 



(2.20) r(i) = 



exp 



E j"i(0<y(0 - ^(W*) + \ j\i{a;T)da 



The equality l|2.19[) follows from this. A straightforward calculation based on H2.20fl , 
(|2~BI) and the Ito formula gives dT(t) = i()'(t)dB(t), where ip(t) is as in l(2~T%j) . Thus 
the proposition follows. □ 

Recall that we have assumed a € (— oo, 1) \ {0} and (| 2 . 1 4|> . Here is the solution 
to Problem PI. 

Theorem 2.3. For T S (0,oo), the strategy (7Pr(i))o<t<T £ -At defined by 

(2.2i) Tr T (t) : = (a')- 1 (t) [(i - /9){A(t) - £(*)} - KW; r)£(<) + i(*M*; r)] 

is i/ie unique optimal strategy for Problem PI. The value function V(T) = V(T, a) 
in (PI) is given by 



(2.22) V(T) = ~[xS (T)] a exp 



(1-a) 



E- =1 / 9j(t;T)dt 



Proof. By Theorem 7.6 in Karatzas and Shreve |21l Chapter 3], the unique optimal 
strategy 7Ty(i) for Problem PI is given by 

ir T (t) := {a'r\t) [^(^(t) + X(t) - £(tj\ , < t < T, 

which, by l|2.18|l . is equal to ^x(i)- Thus the first assertion follows. By the same 
theorem in 21 , V(T) = a- 1 [xS (T)] a r 1 - ct (0). This and {2~T5|> give QFty . □ 

7 



Remark 2.4. We can regard £(t) = J Q * fc(t, s)dY(s), which is the only random term 
on the right-hand side of H2.21JI . as representing the memory effect. To illustrate 
this point, suppose that (otj(t)) is a constant matrix. Then, by (|C.2|) in Appendix 
C, we can express Y(t), whence £(t), in terms of the past prices 5(u), u <E [0,t], of 
the risky assets. 

Remark 2.5. From [211 Theorem 7.6], we also find that 

x ^ T{t)=x so^m 0<t<T 

Remark 2.6. Regarding J23J, we assume this to ensure the existence of solution 
to l|2.10[l for j = 1, . . . , n. Under the weaker assumption l|1.4[) . we could show by a 
different argument that, for j = 1, . . . , n, l|2.10[l has a solution if a 6 (— oo, Oj)\{0}, 
where ay £ (0, 1] is defined by 



(Pj + Qj 



aj := 1 if < Pj < oo, := /2 ( 0) + g 2 if ~ % < Pj < °- 

From this, we see that the same result as Theorem 12 . 31 holds under <|1.4[l if — oo < 
a < a, a ^ 0, where a := minjaj : j = 1, . . . , n}. However, we did not succeed 
in extending the result to the most general case — oo < a < 1, a ^ 0. Such an 
extension, if possible, would lead us to the solution of Problem P3 under 11.4[> (see 
Remark ET%|) . 

3. Optimal investment over the infinite horizon 

In this section, we consider the infinite horizon optimization problem P2 for 
the financial market model M.. Throughout this section, we assume l|2.1[l and the 
following two conditions: 

1 f T 

(3.1) lim - / r(t)dt = f with fe[0,oo), 

(3.2) HmA(t)=A with A = (Ai, . . . , A„)' G R n . 

t — ^oo 

Here recall X(t) = (Ai(t), . . . , A„(i))' from 12. 2[) . In the main result of this section 
( Theorem 13. 40 . we will also assume a* < a < 1, a ^ 0, where 

(3.3) a* := max(a*, . . . , a* ) 
with 

- oo if < pj < 2qj , 



(3.4) 



— 3 — if 2o, < pj < oo. 



Notice that a* e [— oo, -3). 

To give the solution to Problem P2, we take the following steps: 

(i) For the value function V(T) = V(T, a) in (PI), we calculate the following 
limit explicitly: 

(3.5) J(a):= lim -L\og[aV(T)}. 



(ii) For 7r G A in (|3.14|) below, we calculate the growth rate 

(3.6) J*(a) := lim ±- \ogE [(X'>*(T)) a ] , 

I — *oo (X 1 

and verify that J* (a) — J (a). 

(iii) Since the definition of V(T) implies 

(3.7) ]imsup — logE[(X x ' v (T)) a ] < J (a), Vtt G A, 

we conclude that fr is an optimal strategy for Problem P2 and that the 
optimal growth rate J (a) in (P2) is given by J (a) — J* {a) = J (a). 

Let a G (— oo, 1) \ {0} and (3 be its conjugate exponent as in Section [21 For 
j = 1, ... ,n, recall bj(t) from l|2.9[l . We have lim^^ bj(t) = bj, where 

bj ■■= -(l-0)pj -qj. 
Notice that bj < 0. We consider the equation 

(3.8) p)x 2 - 2b jX - (3(1 -/3) = 0. 

When pj — 0, we write Rj for the unique solution /3(1 — (3)/(2qj) of this linear 
equation. If pj > 0, then 

1) + (3(1 - (3) P ) = (i - p)[(pj + Qj) 2 - <$] + q) > q) > o, 

so that we may write Rj for the larger solution to the quadratic equation (|3.8|) . Let 



Kj := ^/b] + (3(l-(3)pl 



Then bj - p]Rj = -Kj < 0. 

As in Sectional we write Rj(t) = Rj(t;T) for the unique solution to p.lOJI in 
the sense of Lemma 12.11 Recall A from H2.17|l . The next proposition provides the 
necessary results on the asymptotic behavior of Rj(t;T). 

Proposition 3.1. Let — oo < a < 1, a / 0, and j G {1, . . . , n} . Then 

(i) Rj(t;T) is bounded in A. 

(ii) liniT-t^oo, t ~>oc Rj(t;T) = Rj. 

(iii) For 5, e G (0, oo) such that S + e < 1, 

lim sup \Rj(t;T) -Rj\ = 0. 

T^oc ST<t<{l-e)T 

Proof. If pj — 0, then lj(t) = and bj(t) — —qj < for t > 0, so that the assertions 
follow from Theorem IB . 31 in Appendix B. 
We assume pj > 0. Since 

the function lj(t) converges to pj exponentially fast as t — > oo. Hence the coefficients 
of the equation (|2.10() converge to their counterparts in (|3.8() exponentially fast, too. 
If — oo < a < 0, then the desired assertions follow from Theorem IB. II in Appendix 
B (due to Nagai and Peng |22j)- Suppose < a < 1. Let ai(t), 02(f) and 03 be as 
in QZD) . Since Rj(t;T) > bj(t)/lj(tf and bj(t)/lj(t) 2 is bounded from below in 
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A, so is Rj(t;T). To show that Rj(t;T) is bounded from above in A, we consider 
the solution Mj(t) = Mj(t; T) to the linear equation 

Mj{t) + 2[a 2 {t) - RjCuityMjit) + a 3 + a x (t)R 2 j = 0, < t < T, Mj(T) = 0. 

Since M 3 (T) - Rj(T) = and 

[Mj(t) - Rj(t)] + 2[oa(*) - Rjai(t)][Mj(t) - Rj(t)] = -oi(t) [Rj(t) - Rj] 2 < 0, 

we have Rj(t;T) < Mj(t;T) in A. However, a 2 (t) - Rjai(t) -» bj - Rjp? < as 
t — > oo, so that Mj(t;T) is bounded from above in A, whence so is Rj(t;T). The 
desired assertions now follow from Theorem lB.2l in Appendix B. □ 

Let j G {1, . . . , n}. For pj{t) in 12.8(1 . we have lim^oo pj(t) — pj, where 



Let Vj(t) = Vj(t;T) be the solution to l|2.14|l as in Section [21 Define Vj by 



Proof. The coefficients of 1(2.1411 converge to their counterparts in 1(3.9(1 . Also, 



Pj := -ppjXj. 



(3.9) (bj - p 2 R 3 ) Vj + 0(1 - P)\ 3 - RjPj = 0. 



Proposition 3.2. Let —oo < a < 1, a / 0, and j e {1, . . . , n} . Then 

(i) Vj(t]T) is bounded in A. 

(ii) lim T ^ t ^ 00i t^oo Vj(t;T) =Vj. 

(iii) For 5, e € (0, oo) such that S + e < 1, 



lim sup \vj(t;T) — Vj\ = 0. 

T^oo ST<t<(l-e)T 



lim [6j(t) - l 2 At)RAt; T)} = bj - p^Rj = -Kj < 0. 

-»oo, t— too J J 



Thus the proposition follows from Theorem IB. 31 in Appendix B. 



□ 



For j = 1, . . . , n and — oo < a < 1, a 0, we put 




Proposition 3.3. Let — oo < a < 1, a / 0. Then the limit J(a) in 1(3.5(1 exists 
and is given by 



(3.12) 




where 



g 3 := v 2 p 2 + 2pjVj - p)Rj -(3(1- (3)~X 2 , j = 1, . . . , n. 



More explicitly, 



(3.13) 
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Proof. Recall gj(t;T) from QTTty . By Propositions IO and EH1 

1 f T 

Jim 7f / 9j{t; T)dt = g jt j = l,...,n. 

T^oo 1 J Q 

From this and J222J, 

-L log [aV(T)] = ^ + i [ rW ^t}J o T T)dt 

n 

1 — a \ _ 

r + -g— " 5j as T -> oo, 

j=i 

which implies l|3.12[l . 

We have % = /3Aj(l - + pjRj)/Kj. Also, 

Pp^PjRj) 2 - 2Pp 2 j R j K j = PfiRjtfjRj - 2Kj) = l/rlijh, - K ) 

= f 3(b 2 -K]) = -/3 2 (l-[3)p 2 j , 

and /3(1 -/3) = -a/ (I - a) 2 . Thus 

v)p) + 2p j v j - /3(1 - p)\ 2 

(3\ 2 

3 

= [{Pp 2 (P3 R s) 2 ~ tPPjRjKj} + 2/8(1 - /3>i(p^i - Ji}) 

+/3p 2 (l - [3) 2 - (1 - (3)K 2 ] 

i(l i,A? [-/3 2 p 2 + 2/3pA- +/3(l-/% 2 ~ {^ 2 + /3(l-/3)p 2 }] 



#f LJ (1-a) 2 K] r 

This and p 2 Rj = bj + ifj imply 



3 = Jl~af K 2 J ~ 
Since (1 — a)(l — f3) = 1, it follows that 

(1 - = (1 - a)[(f3 - l)pj - qj] = -pj + (a - l)q 3 = q 3 a - (p 3 + q 3 ). 

Also, 

K j = (Pj +I3) 2 ~ $V 3 (V3 + 2 %) = (Pi + Qj) 2 + T~—P3(Pj + 2 9j)- 
J 1 — a 

Combining, we obtain (|3 . 1 3|> - □ 

Recall £(t) from 1)2.3(1 . Taking into account H2.21J1 . we consider ir = (ft(t))t>o G -4 
defined by 

(3.14) Ht):=(</)-\t)[(l-0){\(t)-t(t)}-pRZ(t)+p€], t > 0, 

where p := diag(pi, . . . ,p n ) as in Section |2 and R := diag(^i, . . . , R n ), v := 
(ui, . . .,«„)'. 
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Recall that we have assumed ()2.1|) . (|3.1[) and (|3.2|) . Recall also a* from (|3.3|) 
with 1|3.4[1 . Here is the solution to Problem P2. 

Theorem 3.4. Let a* < a < 1. a ^ 0. T/ien 7r is an optimal strategy for Problem 
P2 with limit in (|3.6() . The optimal growth rate J (a) in (P2) is given by 



(3.15) 



J(a) 



^ n ^ n 

2a ^ n ' 2a ^ n ' 



3=1 



3=1 



where Fj 's and Gj 's are as in (|3.10H erne? 13.11|l . respectively. 

Proof. For simplicity, we put X(t) := X x >*(t). For J (a) in lESJ and J* (a) in 
we claim J* (a) = J(a), that is, 



(3.16) 



lim \ogE[X a {T)] = J(a) 
T— too ai 



As mentioned before, (|3.7|l and H3.16|l imply that 7r is an optimizer for Problem P2. 
The equality l|3.15[) follows from this and l|3.13[) 

We complete the proof of the theorem by proving (|3.16l) . 
Step 1. We calculate E[X a (T)]. Define the R-valued martingale L{t) by 



Lit) := exp 



a / {a'(s)w{s)}'dB(s) 



a 



|f7'( S )^(s)|| 2 ds 



t > 0. 



From we have X a (t) = [xS {t)] a L(t) exp[/ Q * N(s)ds] for t > 0, where 



N(t) := a{a'(t)Tt(t)}' 



A(t)-e(t) + -(a-iy(t)7r(*) 



a(l — a) 



W(t)n(t)y [(i - (3){\(t) - am + pRm - h 



= -§ [AW - - (i - -pe}]' 

• [A(t)-$(t) + (l-a){p^(t)-p0}] 

= -~ [{A(i) - £(*)}'{*(<) - ec*)} - (i - «) 2 {pW) - - m] 

Notice that we have used (1 — a)(l — j3) — 1, a(l — j3) = —(3. We write 

1 In 

= -^emm ~ h\mt) + - £. =i Ui (t), 

where 

u i (t):=a(a-l)[^ 2 -(l-/3) 2 A|(t)], i > 0, j = l,...,n, 
and Q = diag(Qx, ...,Q n ), h{t) = h(t; T) = {h^t; T), . . . , hj(t; T))' with 
Qj ■= 13 [1 - (1 - afpjitj] , t>0, j = l,...,n, 
hj{t) := a(a - l)PjRji>j - p\j(t), t>0, j = 1, . . . ,n. 

Therefore, 

E[X a (T)} = [xSo(T)] a ex P '~" 

(3.17) 



1 w /- T 

2 £3=1/0 



x E 



exp 



C(t)Qe(<) + h'(t)at) dt 
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where E denotes the expectation with respect to the probability measure P on 
(O, T T ) such that dP/dP = L(T). 

Step 2. We continue the calculation of E[X a (T)]. We are about to apply Theorem 
[Q in Appendix A to l|3~T7jl . Write B(t) := B(t)-a f* a' (s)ir(s)ds for t > 0. Then 
B(t) is an R n - valued standard Brownian motion under P. By (|2.6[1 . the process 
evolves according to the n-dimensional stochastic differential equation 

(3.18) d£(t) = [y(t) + d(t)£(t)]dt + l(t)dB(t), t > 0, 
where d(t) = diag(rfi(t), . . . , d n (t)), j(t) = diag(7i(t), ...,j n (t)) with 

dj{t) := bj(t) — apjRjlj(t), t > 0, j = 1, . . . , n, 
7j (*) : = Pj (*) + a M (*)% > * > °> J = 1, • • • , n. 

For j = 1, . . . , n, let £/j(i) = Uj(t; T) be the unique solution to the one-dimensional 
backward Riccati equation 

(3.19) U j (t)-l?(t)Uj(t)+2d j (t)U j (t) + Q j =0, < t < T, Uj(T) = 

in the sense of Lemma 13.51 below, and let m,j(t) = rrij(t;T) be the solution to the 
one-dimensional linear equation 

ri»i(t) + [dj(t) - l 2 3 (t)U 3 (t; T)] mj (t) - hj(t) - Uj(t; TfriW = 
^ ' ' < t < T, mj(T) = 0. 

Then, from (|3~l7jl - (j3~27Hl and Theorem lA~Tl we obtain 



(3.21) E[X a (T)} = [xS {T)] a exp 



2 Jo 



T 

fi(t;T)dt 



where, for (t, T) £ A and j = 1, . . . , n, 

fj(t;T) := l^t) m 2 ] (t;T) + 2j J (t)m J (t;T)^l^t)U J (t;T)+u J (t). 
Step 3. We compute the limit J*(a) in l|3.6l) . Let j 6 {1, . . . , n}. Write 

Then dj(t) converges to dj, as t — > oo, exponentially fast. Now 
+PjQj = {bj - ap 2 R 3 ) 2 +p 2 l3 [1 - (1 - a) 2 p 2 i? 2 ] 

= fe 2 - 2aB J (6 j + Kj) + a 2 {b 3 + Kj) 2 +p 2 (3 - a(a - l)(bj + Kjf 
= (1 - a)b 2 + aK 2 + p 2 /3 = b 2 + p 2 (i{\ - 0), 

which implies 



(3.22) s Jd?+plQ j =K j >0. 

Thus we may write Uj for the larger (resp. unique) solution of the following equation 
when pj > (resp. pj = 0): 

(3.23) p 2 x 2 - 2djX - Qj = 0. 

From (|3.22|) . we also see that dj — p 2 Uj = —Kj. Let fhj be the solution to 

(3.24) (dj - p 2 j(Jj) fhj - hj - Ufi, = 0, 
where 

hj := a(a — l)p 2 RjVj — (3\j, *fj := pj + ap 2 Vj. 
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By <|3.21[) . we have 

However, Propositions 13 . 61 and 13 . 71 below imply 

I™ ^ I* f j (t;T)dt = f j 

with 

£ := p 2 rh 2 + 2^m 3 - fiUj + a(a - 1) [p)v 2 - (1 - /3) 2 A 2 ] , 

so that 

1 ™ 

(3.25) J »= f +^E£' 

3=1 

Sfep 4. Here we show that in fact (|3.16(l holds. First, 

pjUj = d + Kj = bj - apjRj + Kj = p?Rj - apjRj = (1 - a)pj Rj , 
whence Uj = (1 — a)-Rj (which we can directly check when p 3 — 0). Next, 
hj + Ujjj = a(a — l)p^RjVj — /3\j + (1 — a)Rj [— 0Pj\j + ccp 2 Vj] 
= Aji-0 + ap 3 R 3 ) = -(1 - a)(3\ 3 [(l - P)+ P jRj}, 

so that 

mj = ^-^/3A,[(1 - f3) + p j R j ] = (1 - a)v 3 . 

K 3 

Therefore, 

h = (1 - afp)v) + 2(1 - a) (ft + afivfivj - (1 - a)p 2 i^ 

+ a(a-l)[^2_ (1 _ /3) 2 X 2] 

= (1 - a) [p 2 v 2 + 2p 3 v 3 - p]R 3 - - 0(1 - /3)A 2 ] = (1 - a)&. 

From (|3~T2l) . l|3~23|l and this, we obtain J* (a) = J(a) or (|3~lU|) . as desired. □ 

In the proof above, we needed the following results. 

Lemma 3.5. Let j S {!,-•■, n}. 

(i) //pj = ; then H3.19(l /ias a unique solution U 3 (t) = Uj(t;T). 

(ii) If Pj > and a* < a < 7 iften <|3.19[l has a unique nonnegative solution 
Uj{t) = UifrT). 

(iii) Ifpj > andO < a < 1, then ()3.19J) /ias a unique solution Uj(t) = U 3 (t;T) 
such that Uj(t;T) > (1 - a)Rj(t;T) forte [0,T], where R 3 (t) = Rj(t;T) 
is the solution to (|2.1()(l in the sense of Lemma XZ.W (iii). 

Proof, (i) When = 0, Ij3.19|l is linear, whence it has a unique solution. 

(ii) For p 3 > and a < 0, we put f(x) = p 2 x 2 — 2b 3 x - /3(1 — /?). Since bj < 
and /3(1 — /3) > 0, the larger solution ^ to /(x) = satisfies p 2 -R 2 < (1 — /3) 2 if and 
only if /((l — (3)/pj) > 0. However, this is equivalent to — 3pj — 2q 3 < (pj — 2q 3 )a. 
Thus, if p 3 > and a* < a < 0, then p 2 R 2 < (1 - /?) 2 or Qj > 0, so that the 
Riccati equation l|3.19|l has a unique nonnegative solution. 
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(iii) Suppose p 3 > and < a < 1. For the solution Rj(t) = R 3 (t;T) to (|2.1U|) 
in the sense of Lemma 12. II (in), we consider 

Let dj(t) be as above. Then, H3.19[l becomes 

Pj(t) - (1 - a)l]{t)P]{t) - 2 [(1 - a^m^t) - dj(t)] Pj(t) 

+ allj^Rjit) -pjRj} 2 = 0, < t < T, 

with Pj(T) = 0. Since (1 - a)l?(t) > and a[lj(t)Rj(t) ~P 3 R 3 } 2 > 0, this Riccati 
equation has a unique nonnegative solution. Thus the assertion follows. □ 

Proposition 3.6. Let «' < a < 1, a / 0, and j 6 {1, . . . ,n}. Let Uj(t;T) be the 
unique solution to (|3.19|) in the sense of Lemma \3.5l and let Uj be the larger (resp. 
unique) solution to (|3.23|1 when pj > (resp. p 3 = 0). Then 

(i) Uj(t;T) is bounded in A. 

(ii) limT^t^oo, t^oo Uj(t;T) = Uj. 

(iii) For S, e G (0, oo) such that S + e < 1, 

lim sup \U 3 (t;T)-U 3 \=Q. 

°° ST<t<(l-e)T 

Proof. We assume < a < 1 and > 0. Since Uj(t;T) > (1 - a)R 3 {t]T) in 
A and R 3 (t;T) is bounded from below by Proposition 13.11 so is U 3 {t;T). Let 
iVj (i) = iVj (i; T) be the solution to the linear equation 

Nj(t) + 2[dj(t) - l?(t)Uj]Nj(t) + Qj + l 2 (t)U 2 = 0, < t < T, Nj(T) = 0. 

By (|3~^2l) . dj(t) - l 2 (t)U 3 -> J, - p?E/j- = -i^- < as t -> oo, so that N 3 {t;T) is 
bounded from above in A. Since Nj(T) — Uj(T) = and 

[Nj(t) - Uj(t)] + 2[dj(t) - l 2 (t)U 3 ][N 3 {t) - U 3 (t)} = -l%t) [U 3 (t) ~ U 3 ] 2 < 0, 

we have, as in the proof of PropositionO U 3 {t; T) < N 3 (t; T) in A. Thus U 3 (t; T) 
is also bounded from above in A. Combining, U 3 (t\T) is bounded in A. The rest 
of the proof is similar to that of Proposition ^. II whence we omit it. □ 

Proposition 3.7. Let a* < a < 1, a / 0, and j G {1, ...,n}. Let m 3 (t;T) and 
ifij be the solutions to (|3.20() and (I3.24H . respectively. Then 

(i) rrij(t;T) is bounded in A. 

(ii) limT-t^oo, t ->oo m-j it; T) = fh 3 . 

(iii) For 5, e £ (0, oo) such that 5 + e < 1, 

lim sup \m 3 (t]T) — rh 3 \ = 0. 

T-»oo ST<t<{l-e)T 

The proof of Proposition 13 . 71 is similar to that of Proposition 13. 21 so we omit it. 

Remark 3.8. We note that the proof of Lemma [3.51 (iii) is still valid under (|1.4J) 
if there were a solution Rj(t) = R 3 (t\T) to H2.10(l . This implies that, to prove an 
analogue of Theorem l3.4l with < a < 1, which is relevant to Problem P3, for 11.4JI . 
one may show the existence of such R 3 (t) when < a < 1. We did not succeed in 
such an extension to Lemma (hi) (see Remark |2.6|) . 

is 



4. Large deviations probability control 

In this section, we study the large deviations probability control problem P3 for 
the market model M. Throughout this section, we assume 1)3. 1[) . (|3.2|) and 

(4.1) either A / (0, 0)' or (pi,..., Pn ) £ (0,...,0). 

For x 6 (0, oo ) and tt G A, let L x ' n (T) be the growth rate defined by 

^(T):= l0S 7 (T) , T>0. 

We have P{L X ^{T) > c) = P (X x ' n (T) > e cT ). Following Pham [231111, we con- 
sider the optimal logarithmic moment generating function 

A(a) := sup limsuplog^exp^TZ^T))], < a < 1. 

irtA T—too 

Since A(a) = aJ(a) for a S (0, 1), it follows from Theorem 13.41 that 

_^ n 1 n 

A(a) =ra+-^,F j {a) + -J2 G,-(a), < a < 1, 
i=i j'=i 

where Fj's and Gj's are as in (|3.10() and (|3 . 1 1 1> . respectively. 

Proposition 4.1. We have (dA/da)(0+) = c and \im a ^i(dA/da)(a) = oo, where 

4 ^ Pi + * 2 

Proof. For < a < 1, -Fj(a) is equal to 

(Pj+gj) 2 * 2 (Pj+gj) 2 A 2 gj 2 a 

[(1 - a){pj + q/) 2 + a Pj ( Pj + 2 qj )} [(l _ a )( Pj + q f + a Pj { Pj + 2 qj )} 2 ' 

From this, Fj(Q+) = A|. This also shows that 

^(a)~A?(l-a)- 2 , afl 
if Pj = and A 3 ^ 0. On the other hand, for < a < 1, 

-^■(a) - -g,- + ^ ^ [(1 - + % ) 2 + apjfa + 2q 3 )] 

^(l-a)Va 
-I 1 

2 [(1 - a)(pj + 9 ,) 2 + apjipj + 2q 3 )] 1/2 ' 
This gives (dGj/da)(0+) = P V[2{ P j + qj)]. This also yields 



aa 2 

if > 0. Thus the proposition follows. □ 
Remark 4.2. From the proof of Proposition ^. II we see that 

dA, . (l-a)- 1 / 2 ^" r- — 

_( a ) „ i J J-.^ ^J Pj (pj + 2 qj ), a T 1 
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if pj > for all j = 1, . . . , n, otherwise 

dA (1-g)- 2 ^ - 2 

^(«) 2— «TL 

For a G (0,1), we denote by 7t(t;a) the optimal strategy in 13.14[1 . Recall 
/(c) from (P3). From Theorem 13.41 Proposition 14. II and Pham 1281 Theorem 3.1], 
we immediately obtain the following solution to Problem P3: 

Theorem 4.3. We have 

1(c) = — sup [ac — A(a)] , c G R. 
ae(o,i) 

Moreover, if a(d) G (0, 1) is such that A(a(d)) = d G (c, oo), then, for c > c, the 
sequence of strategies 

is nearly optimal in the sense that 

lim limsup - logP (x x ^ m (T) > e cT ) = 7(c), c > c. 

Remark 4.4. Theorem 3.1 in Pham is stated for a model different from M. but 
the arguments there are so general that we can prove Theorem l4.3l in the same way. 

We turn to the problem of deriving an optimal strategy, rather than a nearly 
optimal sequence, for the problem (P3) when c < c. We define no G A by 

*„(*) :=(0 _1 («)[A(t)-e(t)], t>0, 
where recall £(t) from l|2.3jl . From l|2.4jl . 

L *,* 0(T)= ^ + 1 f T r(t)dt+ ±. r T i m -0)fdt 



Proposition 4.5. 77ie rate L X ^°(T) converges to c, as T — > oo, m probability. 

Proof. In this proof, we denote by C positive constants, which may not be neces- 
sarily equal. 

For j = 1, . . . , n, we write 

Ai(t) - = [Aj(t) - Xj] + [Xj - Pj K(t)] + N(t), 

where K(t) = J* e^+^-^dB^s) and N(t) = /* e -fe+*)(*- s ) f(s)dB j (s) with 

The process if(i), the dynamics of which are given by 

dK(t) = -( Pj + q,j)K{t)dt + dB 3 (t), 

is a positively recurrent one-dimensional diffusion process with speed measure 
m(dx) — 2e~( pi+q ^ x dx. By the ergodic theorem (cf. Rogers and Williams [311 
v.53]), we have 

Urn - / [Xj - p 3 K(t)] 2 dt = / (Xj - p 3 y)Mdy) = Xj + T( \ , 
T— too 1 J J_ oc Z(p J +q :j ) 
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where v{dy) is the Gaussian measure with mean and variance \/[2{pj + gj)]. 
Since < f(s) < Ce~ 2 ^ 8 , we have 

E [N 2 (t)] < C [ e- 2q ^ t+s Us < Ce- 2q i\ t > 0. 
Jo 

Also, E[K 2 {t)\ < C for t > 0. Therefore, 

\ f E [\{Xj - Pj K(t)}N(t)\] dt <^J^ EiN^t^dt - 0, T —> oo. 
Similarly, 

lim i / [X j (t)-X j ] 2 dt= lim i [ T EliXj^-X^Xj-p.-Kit^dt 
= lim i ( T E [N 2 (i)\ dt = lim I T F, [(A,(i) - X 3 )N(t)} dt = 0. 

T^oo J J T->oo J Jo 

Combining, 

7f / [Aj(t) - &(t)f <ft -> A 2 + ^ T -> oo, in probability. 

Finally, for j = 1, . . . , n and t > 0, 

{*,•(*) ~ 6" W} 2 1 < 2A 2 (t) + 2F [£ 2 (i)] < C 



t 

2/ 



1+ / rJt,s)ds 



I) 



<c, 



so that (1/T) J T [Aj (t) - ^(t^dB^t) -> 0, as T -> oo, in L 2 (ft), whence in proba- 
bility. Thus the proposition follows. □ 

Theorem 4.6. For c < c, ttq is optimal for Problem P3 with limit 

lim - log P (X x ^° (T) > e cT ) = 7(c), c < c. 
Proof. Proposition 14.51 implies limT->oo P (L X ^°(T) > c) = 1 for c < c, so that 
lim - logP (L X <*°(T) > c) = > sup lim - logP (L X '*(T) > c) , c < c. 

T— >oo T Tr£AT—>°° T 

Thus 7To is optimal if c < c. □ 

Remark 4.7. From Theorem 10.1 in Karatzas and Shreve |21l Chapter 3], we see 
that 7To is the log-optimal or growth optimal strategy in the sense that 

sup lim sup - log X x ^ (T) = lim sup - log X x ^ a (T) a.s. 

7rEA T — too 

We also find that lim Q jo 7r(t; a) = 7ro(i) a-S- for t > 0. 

Appendix A. A Cameron-Martin type formula 

In this appendix, we prove a generalization of the Cameron-Martin formula that 
we need in the proofs of Proposition 12.21 and Theorem 13.41 We refer to Myers 
for earlier work. 

Let T € (0, oo) and let R nxn be the set of n x n real matrices. We say that 
A : [0, T] — > R nxn is symmetric if A(t) is a symmetric matrix for all t G [0, T]. Let 
(ft, T, P) be the underlying complete probability space equipped with filtration 
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(•7"t)o<t<T satisfying the usual conditions. We assume that the R"-valued process 
£(t) satisfies the n-dimensional stochastic differential equation 

(f) = [ a (t) + b(t)£(t)]dt + c(t)dB(t), < t < T, 

where B(t) is an R™- valued standard (^ r t )-Brownian motion and all the coefficients 
a : [0, T] — > R™ and b, c : [0, T] — > R nxn are deterministic, bounded measurable 
functions. 

Theorem A.l. Let Q : [0,T] -> R nx " and h : [0,T] -> R n 6e deterministic, 
bounded measurable functions. We assume that Q is symmetric. We also assume 
that there exists a bounded symmetric function R : [0,T] — > R nxn satisfying the 
backward matrix Riccati equation 

R(t) - R(t)c(t)c'(t)R(t) + b'{t)R{t) + R(t)b(t) + Q(t) = 0, < t < T, 
R(T) = 0. 

Let v : [0, T] —> R™ be the solution to the linear equation 

v(t) + [b(t) - c{t)c'{t)R{t)]'v{t) - h(t) - R(t)a(t) = 0, < t < T, 
v(T) = 0. 
Then, forte [0,T], 



(A.l) 



(A.2) 



E 



exp ■ 



'-X'{s)Q{s)i{s)+h'{s)i{s))ds 



= exp 



v'(tm ^'(t)Rmt) 



x exp 



{v'{s)c(s)c'(s)v(s) + 2a'(s)v(s) - tr (c(s) c' (s)R(s))}ds 



Proof. We put K(t) = c'(t)[v(t) - R(t)£(t)] for t e [0, T]. Then, by the ltd formula, 



d 



v'(m) - \e(t)R(t)m 



K'{t)dB{t)--\\K{t)fdt 



-^'(t){R{t) + R(t)b{t) + b'{t)R{t) - R{t)c(t)c' {t)R(t)}i{t) 
+ {v{t) + (b(t) - c{t)c\t)R{t))'v(t) - R(t)a{t)}'i(t) 

+ ^{v'(t)c(t)c'(t)v(t) + 2a'{t)v(t) - tv(c{t)c'{t)R{t))} 



+ -{v'(t)c(t)c'(t)v(t) + 2a'(t)v(t) - tr(c(t)c'(t)R(t))} 
Therefore, K'{s)dB s — | \\K(s)\\ 2 ds is equal to 



dt 



dt. 



1 f T 

- o / W(s)c(s)c'(s)v(s) + 2a'(s)v(s) - tv(c(s)c'(s)R(s))} ds. 

2 Jt 
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Since K(t) is a continuous Gaussian process, the process 

M{t) := cxp Ij K'(s)dB(s) ~\ j \\K{s)fds^ , < t <T, 
is a martingale (cf. Example 3(a) in 23J Section 6.2]). Thus 

explj* K\s)dB{s)-\£ \\K(s)\\ 2 ds ] 



E 



Combining, we obtain the theorem. □ 

Appendix B. Asymptotics for a solution to Riccati equation 

Here we summarize the results on the asymptotics for a solution to Riccati or 
linear equation that we need in Section [3J 

For T G (0, oo), we consider the one-dimensional backward Riccati equation 

(B.l) R{t)- ai (t)R 2 (t) + 2a 2 {t)R{t)+a 3 (t) = 0, < * < T, R(T) = 0, 

where 

(B.2) 0i (.) G C([0, oo) -> R) for i = 1,2, 3, 

(B.3) ai(t) > for f > 0, 

(B.4) for i = 1. 2, 3, a,(t) converges to 3j exponentially fast as i — > oo, 

(B.5) Si > and d 2 2 + axa 3 > 0. 

By (IB. 5(1 . we may write R for the larger solution to the quadratic equation 

axR 2 - 2a 2 R -03 = 0. 

Recall A from (|2~TTll . 

Theorem B.l (Nagai and Peng |27], Section 5). We further assume 
(B.6) a 3 (t) > for t > 0. 

Then, for T G (0, oo), IjB.ip /ias a unique nonnegative solution R(t) = R(t; T), and 
it satisfies the following: 

(i) R(t;T) is bounded in A. 

(ii) limT^t^oo^ i^oo T) = R. 

(hi) For S, e G (0, oo) smc/i t/iat <5 + e < 1, 

lim sup \R(t; T) - R\ = 0. 

T—oo 5T<t<(l-e)T 

When the condition 1B.6|) is lacking, we have the following: 

Theorem B.2. We assume that, for every T G (0, oo), the equation (|B.1|I has a 
solution R(t) = R(t;T) that is bounded in A. Then (i)— (iii) in Tfeeorem lB.il hold. 

The proof of Theorem IB. 21 is almost the same as that of Theorem IB. II in [27] , 
whence we omit it. 

We turn to the one-dimensional backward linear differential equation 

(B.7) v(t) - bi(t;T)v(t) + b 2 (t;T) = 0, < t < T, v(T) = 0, 
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where 

(B.8) k(-;T) G C([0,T] -> R) for T G (0,oo) and i = 1,2, 

(B.9) lim bi(t;T) = bi for i = l,2, 

(B.10) 61 > 0. 

Theorem B.3 (|27|, Section 5). For T G (0, 00), write v(t) = v(t;T) for the 
solution to (|B.7|1 . Let v be the solution of the linear equation b\V — 62 = 0. Then 

(i) v(t;T) is bounded in A. 

(ii) limr-t-foo, t-^ooV^T) =v. 

(hi) .For o", e G (0, 00) smc/i t/ia£ (5 + e < 1, 

lim sup \v(t; T) - v\ = 0. 

T^oo 5T<t<(l-e)T 

Appendix C. Parameter estimation 

In this appendix, we use the special case of our model M. in which <jjj-(i)'s are 
constants, i.e., 

°~ij(t) = o-ij, t>0, i,j =l,...,n. 
We explain how we can statistically estimate the parameters <7y , pi and qi from 
stock price data. This problem, for the univariate case n = 1, is discussed in |3ll20|. 
Here we are interested in the multivariate case n > 2. As for the expected rates 
of return /Zj, there is as usual a structural difficulty in the statistical estimation of 
them (cf. Luenberger [2U Chapter 8]), whence we do not discuss it here. 
From l|1.5|) . we see that 

(C.l) E[Y J 2 (t)]/t = f(t;p j ,q j ), <>0, j = l,...,n, 

where 

/( * ;P ' 9) - = ¥T^ + 1^F * ' >0 

(cf. [I], Examples 4.3 and 4.5). Notice that f(t;0,q) = 1. From or (JlTJ and 
the Ito formula, the solution S(t) = (Si(t), . . . , S n (t))' to (f 1 . 1|> is given by 



(C.2) Si(t) = Siexp 



t > 0, 



for i = 1, . . . , n. Since Y(t) has stationary increments, we may define 

^.(t- s): = J_cov{logf^,log|^}, i> S >0, i,j = l,...,n, 
t-s { Si(s) Sj(s)) 

where cov(-, •) denotes the covariance with respect to the physical probability mea- 
sure P. By l|C.lj) and !|C.2j) . we see that 

En 
&i m o-j m f(t;p m ,q m ), t>0, i,j = l,...,n. 

Suppose that we are given data consisting of closing prices of n assets observed at 
a time interval of N consecutive trading days. For m = l,...,N and i = 1, . . . , n, 
we denote by Si(m) the price of the ith asset on the mth day. Notice that here the 
time unit is the day. Pick M < N, and, for t G {1, . . . , M}, define Ui(m) = Ui(m, t) 

by 

S-' { Tfl ~\~ t ) 

Ui(m) := log , m — 1,2,..., TV t. 

Si(m) 
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For t = 1, 



(C.3) 



, M, we consider the estimator 
252 ^JV-t 



u y (t) := 100 



EiV —t 
{«i(m) -tZi}{«j(m) 



of (t) , where Mj 



f(iV-t-l) 

(TV — t) -1 X)m=i u i{ m )- The number 252, which is the average 



number of trading days in one year, converts the return into that per annum, while 
the number 100 gives the return in percentage. 

We estimate the values of the parameters <iij , pi and qt by nonlinear least squares. 
More precisely, we search for the values of them such that the following least squares 
error is minimized: 

M n n 

t—i i=i j=i 

We show numerical results obtained from the following daily stock prices from 
September 18, 1995, through September 16, 2005: 

Si : Pfizer Inc., S 2 : Wal-Mart Stores Inc., S 3 : Exxon Mobil Corp. 

Here we use closing prices adjusted for dividends and splits, which are available at 
Yahoo! Finance [HSj) rather than actually observed closing prices. In this example, 
we have 

n = 3, N= 2519, M = 100. 
The estimated values of the parameters are as follows: 

on cri2 CT13] [28.7 -14.1 9.1] fpi q{] [0.086 0.305" 
o- 2 i cr 22 cr 23 = 20.4 22.3 13.4 , p 2 q 2 = 0.261 0.044 
C31 cr 3 2 cr 33 -1.8 - 4.6 24.9 p 3 q 3 0.076 0.098 



Using the signed square root SSR[x] := sign(a;)y / [xf, we write 

Dytt) :=SSR[%(t)], dy(t):=SSR[t;y(t)], t > 0, i, j = 1, . . . ,n. 

In Figures IC~T1^. 31 the dotted lines are the graphs of djj(t)'s, while the correspond- 
ing solid lines represent those of Dij(t)'s that are obtained by using the nonlinear 
least squares above. We see that the fitted functions Dij(t) simultaneously ap- 
proximate the corresponding sample values dij(t) well for this data set. We have 
repeated this procedure for various data sets and obtained reasonably good fits in 
most cases. 
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